Effects of Zeeman field on a Spin Bose-Metal phase 
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We consider Zeeman field effects on a Spin Bose-Metal (SBM) phase on a two-leg triangular ladder. This 
phase was found in a spin- 1/2 model with ring exchanges [D. N. Sheng et. al., Phys. Rev. B 79, 205112(2009)], 
and was also proposed to appear in an interacting electronic model with longer-ranged repulsion [Lai et. al., 
Phys. Rev. B 81, 045105 (2010)]. Using bosonization of a spinon-gauge theory, we study the stability of the 
SBM phase and its properties under the field. We also explore phases arising from potential instabilities of the 
SBM; in all cases, we find a gap to spin-1 excitations while spin-nematic correlations are power law. We discuss 
two-dimensional analogues of these phases where spinons can pair with their own species. 

PACS numbers: 



I. INTRODUCTION 

There has been much recent interest in gapless spin liq- 
uids stimulated by several experimental candidates, includ- 
ing two-dimensional (2D) triangular lattice based organic 
compounds '--'^ k-(ET)2Cu2(CN)3 and EtMe3Sb[Pd(dmit)2]2 
and 3D hyper-kagome material^ Na4lr308. One line of theo- 
retical ideas considers states with a Fermi surface of fermionic 
spinonsii^ For the 2D spin liquids, such a state arises as a 
good variational wavefunction^ for a spin model with ring ex- 
changes and is also an appealing candidate for an electronic 
Hubbard model near the Mott transition] ^'"i'^ 

Driven by the need for a controlled theoretical access to 
such phases, Ref. |T3| considered the Heisenberg plus ring ex- 
changes mode l on a two-leg triangular strip - so-called 
zigzag chain. Using numerical and analytical approaches, 
Ref. [Tsl found a ladder descendant of the 2D spin liquid 
in a broad range of parameters and called this phase "Spin 
Bose-Metal" (SBM). The name refers to metal-like itinerancy 
present in the spin degrees of freedom (note that there is no 
electric transport to speak of in the spin-only model). Further 
work Ref. [la studied electronic Hubbard-type models with 
longer-ranged repulsion and showed that they are promising 
systems to realize such an SBM phase in a Mott insulator 
of electrons proximate to a two-band metallic phase on the 
zigzag chain. 

This paper continues efforts to gain insights about the 2D 
spin liquid from the solvable 2-leg ladder example. Here we 
study the SBM phase under Zeeman magnetic field, while 
in a separate paper we will study orbital field. One motiva- 
tion comes from experiments on the 2D spin liquid materials 
K-(ET)2Cu2(CN)3 and EtMe3Sb[Pd(dmit)2]2 measuring ther- 
modynamic, transport, and local magnetic properties under 
strong fieldS] ''-^'^i'^~ — An important question is whether the 
field can induce changes in the physical state of the system. 

To this end, we explore possible instabilities of the 2-leg 
SBM state in the Zeeman field. There have been many stud- 
ies of 2D and ID spin models under magnetic field show- 
ing rich behaviors. For example, the phase diagram of the 
Ji — J2 antiferromagnetic chain with Ji, J2 > in the 
field22r— contains one-component and two-component Lut- 
tinger liquids, a plateau, a phase with static chirality order, 
and a phase with spin-nematic correlations. In the spirit of 



such studies, we allow a large range of fields, which could 
be numerically explored in spin or electronic models realizing 
the SBM phaseji^ii^ We remark that experiments on the spin 
liquid materials achieve only relatively small fields - e.g., the 
maximum magnetization is < O.OI^b per spin. Nevertheless, 
some of our 2-leg ladder phases from the broader theoretical 
study motivate interesting 2D states that are worth exploring. 

The SBM phase on the zigzag chain can be viewed as a 
Gutzwiller-projected spinon state where both '[ and I spinon 
species populate two Fermi segments (cf. Fig.[rii. The projec- 
tion eliminates the overall charge mode leaving three gapless 
modes. We find that this phase can in principle remain stable 
under the Zeeman field. We also identify all possible instabil- 
ities out of the SBM. 

Loosely speaking, the instabilities correspond to pairing of 
spinons separately within each species (a kind of triplet pair- 
ing). More precisely, the relevant interactions can be inter- 
preted as moving a "Cooper pair" from one band to the other 
of the same species. Of course, there is no long-range pairing 
order in the quasi- ID and in fact the dominant correlations 
in our system need not be of "pair-type" - the Bosonization 
provides the proper treatment, while this language is only for 
convenience. 

It can happen that the pairing is relevant for one spinon 
species but not the other In this case the system retains 
two gapless modes. Interestingly, spin-1 excitations become 
gapped (i.e., transverse spin coiTelations are short-ranged), 
while spin-2 excitations are gapless (i.e., nematic or two- 
magnon correlation functions show power law). 

It can also happen that the pairing is relevant for each 
spinon species. In this case the system retains only one gap- 
less mode. Again, spin-1 excitations are gapped while spin-2 
remain gapless. It further turns out that the system breaks 
translational symmetry and has either period 2 Valence Bond 
Solid (VBS) or period 2 static chirality order 

Such thinking about pairing within the same spinon species 
can be extended to 2D. Here, if we pair only one species and 
not the other, we have a gap to spin-1 excitations while at the 
same time we have critical spin-2 correlations and the system 
retains the gapless Fermi surface for the unpaired species. On 
the other hand, if we have pairing within both spinon species, 
the system acquires a long-range spin-nematic order.— 

Spin-nematic phases were discovered and much discussed 
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recently in other interesting frustrated systems. For instance, 
such phases were found in the antiferromagnetic zigzag ladder 
with easy-plane anisotropy2^ and in the ferro/antiferro zigzag 
ladder ( Ji < 0, J2 > 0) in the Zeeman field.-"" As for 
examples in 2D, spin-nematic order was found in the frus- 
trated square lattice with ferromagnetic Ji < and antiferro- 
magnetic J2 > and ring exchanges?^ and in the triangular 
lattice with ferromagnetic Heisenberg and antiferromagnetic 
ring exchanges.i'? Though, many details of the nematic phases 
proximate to the SBM studied here are of course different. 

The paper is organized as follows. In Sec. [Ill we consider an 
electronic Hubbard-type model with longer-ranged repulsion 
under Zeeman magnetic field and discuss the weak coupling 
phase diagram in the two-band regime. We then take a leap 
to the Mott insulator regime, which can be achieved from the 
electronic perspective by gapping out the overall charge mode 
using an eight-fermion Umklapp interaction. In Secs. Hlllirvl 
we discuss the theory and properties of the SBM under Zee- 
man field, and in Sec.[V]we consider possible instabilities and 
characterize the resulting phases. We conclude by discussing 
generalizations of these phases to 2D. 



II. ELECTRONS ON A TWO-LEG ZIGZAG STRIP IN A 
ZEEMAN FIELD: WEAK COUPLING APPROACH 

In this section, we consider half-filled electronic ti — ^2 
chain with extended repulsive interaction in the magnetic Zee- 
man field. The Hamiltonian is H ^ Hq + Hz + Hy with 

Hq = -Y^ [tici{x)Cc,{x + 1) + t2cUx)Ca{x + 2) 

a:, a 

+H.C.] , (1) 

Hz = -hY,S'{x) , (2) 

X 

Hv = ^J2^^^~ x'Mx)n{x') . (3) 

Here Cq {x) is a fermion annihilation operator, x is a site la- 
bel on the one-dimensional (ID) chain, and a =^^, | is a spin 
index; n{x) = c^^{x)c^{x) + c|(x)c^(.t) is electron number 
on the site. Throughout, electrons are at half-filling. The Zee- 
man field couples to electron spin S^{x) = ^[c!^{x)c^{x) — 
c|(x)c;(a;)]. 

In the weak coupling approach, we assume Hy <C Hq, Hz 
and start with the non-interacting band structure given by 
^^0 + Hz and illustrated in Fig. [T] In this paper, we focus 
on the regime t2/ii > 0.5 and not too large Zeeman field so 
that there are two occupied Fermi segments ("bands") for each 
spin species. The corresponding phase boundary in the tijti- 
h/ti plane is shown in Fig.|2] For fields exceeding some crit- 
ical values, the second spin-^ Fermi segment gets completely 
depopulated; this regime leads to a different theory and is not 
considered here. 

The spectrum is linearized near the Fermi points and the 




FIG. 1 ; Single-particle spectrum in the presence of the Zeeman field, 
^^/^(fc) = —2t\ cos(fc) — 2t2 cos(2fc) ^ ^ — II, shown for param- 
eters t2/t-i = 1 and h/t\ = 1/2. Our kp-s. denote right-moving 
momenta G (— t!", t!"); with this convention, the half-filling condition 
reads kpi-x + kpii + kF2-\ + kF2i — — tt. 
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FIG. 2: Free electron phase diagram in the t2/t\-h/t\ plane. In this 
paper, we focus solely on the lower region where both spin species 
have two Fermi seas ("bands"). For reference, we give the magneti- 
zation = {rif — ni) I (n-f + nj,) at the transition for several band 
parameters: M^Ht = 0.32, 0.46, 0.54 for f2/ti = 1.0, 1.5, 2.0. 



electron operators are expanded in terms of continuum fields, 

c„(x)=^e'^'=— cp,„, (4) 

with P = R/L = +/— denoting the right/left movers and 
a = 1 , 2 denoting the two Fermi seas for each spin species, 
cf. Fig.[T] There are four different Fermi velocities Vaa- 

Using symmetry arguments, we can write down the most 
general form of the four-fermion interactions which mix the 
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right and left moving fields: 

-ffint = H-f + H_i+ H-fi , (5) 

iPRlaPLla + ^22PR2aPL2a (6) 

+ ^12{PR.laPL2a + PLlaPR2a) O) 

+ w"2(4laCllaCL2aCi?2a +H.C.) , (8) 

= ^^li{PRatPLbl + PLatPRbl) ■ (9) 



(Interactions that do not mix right and left movers only shift 
velocities and do not affect the weak coupling treatment.) 
The weak coupling renormalization group (RG) equations 
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(10) 

(11) 
(12) 

(13) 
(14) 



Here O = dO/dl, where £ is logarithm of the length scale; 
a =t, I; and a, & G {1, 2}. We see that the terms A^j; do not 
flow and the two spin species behave independently from each 
other in the weak coupling regime. 

We therefore focus on one species at a time. Effectively, 
this is equivalent to a two-band model of spinless fermions in 
one dimension o'^'-^-^" — in the absence of any Umklapps. The 
RG Eqs. (fT0]i-(fT3Tl have the Kosterlitz-Thouless form and can 
be solved exactly. We define 
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^11 _^ 

27r?;iQ 27rU2a T^ivia + V2a) 



Eqs. (fT0li-(fT3]l simplify, 

•a _ (via + V2a? + '^VlaV2a ,^„a ,2 
2n''viaV2a[Vla + V2aY 

Wi2 = -y°'w^2 ■ 



(15) 



(16) 
(17) 



The renormalizes to zero if the bare couplings satisfy 



y-ii = 0) > /|^l^±i^!±i!^ = 0)1 .(18) 

In this case, the two-band state of species a is stable and gives 
two gapless modes. 

On the other hand, if the condition Eq. ( fTST i is not satisfied, 
then runs to strong coupling. In this case, only one gap- 
less mode remains. To analyze this, we bosonize 



CPa 



Va 



with canonically conjugate boson fields: 

[ipaaix),ipbi3ix')] = [6*00 (cc) , 6*6/3 (a;')] = , (20) 

[(Paaix),9bi^{x')] = inSabSalsOix - x') , (21) 

where 8(.x) is the Heaviside step function. Here we use Ma- 
jorana fermions {{'i]aa,i]b0} — 25abSap) as Klein factors, 
which assure that the fermion fields with different flavors anti- 
commute with one another 
For convenience, we introduce 
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a =t or I , 




71-f 



72t 
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~ V2 
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(22) 
(23) 
(24) 



and similarly for (p variables. The term becomes 

Wr2{cRlaCLla(^L2aCR2a +ii-C.) ^ W?2 COs(2\/2l^;; ) . (25) 

When is relevant and flows to large values, it pins the 
difference field cp^, while the overall field remains gap- 
less (as it should, since the a-electrons have an incommen- 
surate conserved density and there are no four-fermion Umk- 
lapps). In this phase, the a-electron operator becomes gapped. 
Pair-a-electron operator is gapless, and also specific particle- 
hole composites are gapless, with details depending on the 
sign of wf2- We are primarily interested in repulsively in- 
teracting electrons and expect the particle-hole observables to 
be more prominent, although not dramatically since for too 
strong repulsion the conducting state of the t and I electrons 
is destroyed towards Mott insulator as described below. We 
do not provide more detailed characterization of the conduct- 
ing phases of electrons here, as we are eventually interested in 
the Mott insulating regime where the ^ and I. species become 
strongly coupled. (The two-ba nd sp inless electron system was 
considered, e.g., in Refs. and our analysis in Sec.lIVI 

can be readily tailored to the electronic phases here.) 

In the model with longer-ranged density-density repulsion, 
Eq. ([3]l, the bare couplings are 



All = ^Q=0 - ^2feFi„ , 
A22 yQ=0 - V2fcj^2„ I 

A12 = yQ=o - Vkp^^+kf 



Alt - ^Q=o 



(26) 
(27) 
(28) 
(29) 
(30) 



Here Vq = Y.7'=-oo ^(^ " x')e'^^^-^ > = V^q. 
As an example, we consider the following potential 



V{x-x') 



U 



\x- 
\x- 



= 
> 1 



(31) 



(19) 



This was used in Ref. to provide stable realizations of the 
C2S2 metal and the SBM Mott insulator of electrons in zero 
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FIG. 3: An example of the weak-coupling phase diagram in the 
electron system under the Zeeman field, using model interactions 
Eq. dsn with k — 0.5 and 7 = 0.3. We focus on the region where 
the kinetic energy gives four modes (cf. Figs. [T] O and find four 
phases: metallic phase with four gapless modes evolving out of the 
C2S2 phase in zero field; phase with three gapless modes where only 
the lulj-term is relevant and flows to strong coupling; phase with 
three gapless modes where only the tii^j'te™ is relevant; and phase 
with two gapless modes where both the ^"d ars relevant. 
The uijj-term is relevant in the region with hash lines at roughly 45 
degrees and the uifj-term is relevant in the region with hash lines at 
135 degrees with respect to the horizontal axis. Note that the uijj- 
term always becomes relevant upon approaching the boundary of the 
two-band structure.— 



field. Here U is the overall energy scale and also the on- 
site repulsion; dimensionless parameter k controls the rela- 
tive strength of further-neighbor interactions; and 7 defines 
the decay rate. Applying the stability condition, Eq. ( fTsT i, we 
can now determine the phase diagram in the weak coupling 
approach in the regime where the kinetic energy gives four 
modes. 

Figure |3] provides an illustration for 7 = 0.3 and k = 0.5. 
The interaction is relevant in the region with hash lines 
at roughly 45 degrees and the is relevant in the region 
with hash lines at 135 degrees with respect to the horizontal 
axis. There are four distinct phases. First, when both and 
irrelevant, we have a phase with four gapless modes, 
which is connected to the C2S2 phase at ft, = 0. (Note, how- 
ever, that we assumed Hz ^ Hy, so the formal ft. — ^ limit 
here is different from the weak coupling analysis at ft, = in 
Refs.fllsel) 

Next, when is relevant while u)^2 is irrelevant, we have 
a phase with three gapless modes: one associated with the 
^-electrons and two associated with the ^-electrons. In this 
phase, inserting a single ^-electron costs a finite gap while in- 
serting a pair of ^-electrons or a particle-hole combination of 
^-electrons is gapless. The | electrons are completely gapless. 

When w\2 is relevant while w\2 is irrelevant, we have an- 
other phase with three gapless modes, which is similar to the 
preceding paragraph but with ^ and I interchanged. As can 



be seen in Fig. [3] W12 is always relevant when ft approaches 
the critical valuer^ and the instability arises because the V21 
approaches zero. 

Finally, for large ti, both w\2 and 111^2 are relevant and 
we have a phase with only two gapless modes: one associ- 
ated with spin-f and the other with spin-| species. In this 
case, inserting a single electron of either spin is gapped, while 
inserting a pair or a particle-hole combination of same-spin 
electrons is gapless. 

III. TRANSITION TO MOTT INSULATOR: SBM PHASE 

Note that all phases accessed from the weak coupling anal- 
ysis are conducting along the zigzag chain. Mott insulat- 
ing states do not appear since there is no four-fermion Umk- 
lapp. The half-filled system does become insulating for suf- 
ficiently strong repulsion. This can be achieved by including 
a valid eight-fermion Umklapp, which is irrelevant at weak 
coupling but can become relevant at intermediate to strong 
couphng>i^ii^ 

Hs = U8(4i.,.4^42t42;CLltCLUCL2tCL2; +H.C.) 

- 2t;8cos(40p+) , (32) 

where ^p-^- is defined in Eq. ( |23] l and describes slowly vary- 
ing electron density, Pe{x) = 2dx0p+/iT. The density- 
density repulsion gives coarse-grained interaction Hiat ^ 
VQ=o{dxOp+)'^ . This will stiffen the 9p+ field and will reduce 
the scaling dimension of the Umklapp term. For sufficiently 
strong repulsion the Umklapp becomes relevant and will grow 
at long scales, pinning the 9p+ and driving a metal-insulator 
transition. As discussed in Refs. [l3lfT6i we expect that Mott 
insulator corresponding to a spin model with spins residing on 
sites is described by wg > and the pinning condition 

46'^'^ = TT (mod 27r) . (33) 

Such gapping out of the overall charge mode can occur 
out of any of the four conducting phases discussed in Fig. |3] 
When this happens out of the four-mode metal, we obtain 
spin liquid Mott insulator with three gapless modes - the Spin 
Bose-Metal. In principle, one could perform an intermedi- 
ate coupling analysis similar to that in Ref. [l^ to estimate the 
strength of the repulsion needed to drive the metal-insulator 
transition, but we will not try this here. Below we discuss 
qualitatively the stability and physical observables in the SBM 
phase under the Zeeman field. We will then consider instabil- 
ities of the SBM similar to the ti;"2"driven transitions out of 
the four-mode metal above, but now with the '[ and l systems 
strongly coupjed to form the Mott insulator. 

Reference |j_3| also presented another route to describe the 
SBM in a spin-only model by using Bosonization to analyze 
slave particle gauge theory. The formalism is similar to the 
electron model analysis, but with electron operators Cq(.t) re- 
placed with spinon operators (x) and the gauge theory con- 
straint realized via an explicit mass term for 9p^, 

-^gaugc theory = m (^p+ — ) ■ (34) 
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Loosely speaking, spinons are electrons that shed their overall 
charge once the Umklapp term iJg became relevant. Noter^ 
however, that in the spin-only model, there are no free 
spinons, unlike the situation in the electronic model where we 
have electron excitations above the charge gap. 

From now on, we will use the spinon-gauge language. To 
get some quantitative example, we consider the case where 
spinons do not have any interactions other than Eq. ( l34b . i.e., 
all residual interactions like Eq. Q are set to zero. Once the 
9p^ field is pinned and after integrating out the ipp^ , we obtain 
an effective action for the remaining fields [6^^ , 9^ , 9^ ) = 

&^ and {(fa+, , f^) = defined in Eqs. 



1 

2^ 



Matrix elements of A and B are. 



(35) 
(36) 
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/ ±1 



B 



V2v 



V2v 
2v 



2v 



V2v 
2v 



where 



± Via ±V2a . . 

i-'a = t; 1 a=T,l; 



(37) 



Having all the matrix elements, we can numerically calculate 
the scaling dimensions of the lu^j-terms in Eq. dZSl ). 

As an illustration. Fig. |4] shows the results along a vertical 
cut at t2/ti ~ 1 from Fig.[3](assumed driven into the Mott in- 
sulator as described above). We see that in the absence of the 
residual interactions the SBM remains stable under the Zee- 
man field. We also note that the scaling dimensions of the 
and 11)^2 have opposite trends, which implies that the overall 
stability is reduced. Since the scaling dimension of the 
interaction decreases with increasing field, it is likely that this 
will be the first instability channel upon including the resid- 
ual interactions. This finding is similar to the weak coupling 
analysis where the ^-system tends to become unstable first. 
We want to emphasize, however, that neglecting the residual 
spinon interactions is likely a poor approximation for any real- 
istic spin model, and any calculations in this scheme should be 
taken with caution. The only precise statement here is that the 
SBM can in principle remain stable under the Zeeman field. 

In Sec. [V] we discuss phases proximate to the SBM. Mo- 
tivated by the above observations, we will consider first the 
case where only the term becomes relevant; we will also 
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FIG. 4: Scaling dimensions Aftu^j] ™d A[to|2] as a function of 
h/ti for fixed t2/ti — 1, calculated in the absence of residual spinon 
interactions. In this case, the scaling dimensions stay greater than 2 
and the SBM phase remains stable under the Zeeman field. 



consider the situation where both w{2 and w]^2 relevant. 
Before this, we need to describe main physical observables 
in the SBM under the Zeeman field, which we will then use 
to analyze the instabilities and the properties of the resulting 
phases. 



IV. OBSERVABLES IN THE SBM IN ZEEMAN FIELD 

In the presence of the Zeeman field, the system has S^^^ 
spin conservation symmetry and complex conjugation sym- 
metry (C : i — > — i) in the basis. The system also has lat- 
tice translation and inversion (/ : x — > —x) symmetries. The 
internal symmetries are suffici ently reduced compared with 
the SU(2)-invariant case of Ref. [Ij that we need to revisit the 
physical observables in the SBM. 

We first consider 5*^ -conserving bilinears, which we will 
also call "spin-0" objects. 



— /laa/fliQ J (39) 



1 



/l 1 a 2 a /l 2 a /-R 1 a 



(41) 



2 

2 {^fLlcjL2a + fmafniaj , (42) 
[fllafL2a - fR2afRla ) , (43) 



(no summation over a or a). We define e_Q ~ Eq and x-Q = 

Xq so that e{x) and x(a;) are Hermitian operators. 

The e bilinears appear, e.g., when expressing spinon hop- 
ping energies, while the x bilinears appear in currents. Specif- 
ically, consider a bond [x, x + n], 

- fl{^)Ux + n)+ H.c. , (44) 
^ i [flix)Ux + n) - H.c] , (45) 
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where a ='[ or I species can come with independent ampli- 
tudes. Expansion in terms of the continuum fields gives, up to 
real factors. 



and also 



(n) 

Q 

(») 



(46) 
(47) 



Note that we can view e{x) as a site-centered energy oper- 
ator, e.g., e{x) - 6(1) (.T - 1) + 6(1) (x) - B^^^x - 1), in 
the sense of having the same symmetry properties. We can 
also view e{x) ^ S^{x) in the same sense because of the 
presence of the Zeeman energy. [More generally, the symme- 
try properties of any operator are not changed upon multiply- 
ing by S'^{x).] On the other hand, the bond operator 6("'(x) 
has the same symmetry properties as a bond energy such as 
S{x) ■ S{x + n) and can be used to characterize VBS correla- 
tions in the spin system. 

Similarly, we can view xi^) a site-centered current, 
X{x) ~ J(i)(a; - 1) + j'^'^Hx) - j'-^^x - 1), and also 
as a scalar chirality, ^ 1) • S'(.t) x 5(a; + 1), while 

J^'^\x) has the same symmetry properties as a spin current, 
J^(")(a;) - i[S+{x)S-{x + n)-K.c.]. 

Symmetry analysis shows that eg transforms to e_Q under 
either lattice inversion / or complex conjugation C, while XQ 
transforms to — X-Q under either / or C. We can then give 
an independent argument for the relations Eqs. ( |46l ) and ( |47] i 
for Q 7^ 0, TT, and can show generally that, up to complex 
phase factors, such eg and XQ cover all independent spin-0 
observables for the system in the Zeeman field. 

Special care is needed for Q = it. In this case, Eqs. ( |46] l 
and (l47b hold only for n ~ even. On the other hand, g("-°'i'i) 
is odd under inversion / and even under complex conjugation 
C, while 

is even under / and odd under C. In partic- 
ular, the nearest-neighbor bond i?!^'' and JT^^'' are independent 

(2) (2) 

observables from e^r 6^ and Xtt ~ JTir . In the present 
SBM problem, such Q = tt observables do not appear as bi- 
linears but appear as four-fermion terms below. 

The bosonized expressions for the spin-0 bilinears are: 



=*(''''+ +"^-+)sin(V2^^) ,(49) 
(50) 

-imcmo^e'^'- , (51) 

Vic.mc.e'^'" cos(^/2^-) , (52) 

where we used definitions Eqs. (I22]i-(l24li and a = +/— for 
spin t or 4 and a = +/— for band 1 or 2. 

To bring out the wavevector Q — ir that will play an impor- 
tant role in the analysis of phases near the SBM, we need to 
consider four-fermion terms. We find, 

6^1): i(efcpi^+fc^2t^fcFu+fcF24, ~ H.c.) (53) 
- fsin(\/2(/j::^)sin(V2(p^)sin(2e'p+); (54) 

«(XfeFit+fcF2TXfeFu+fcF2i - ^ (55) 
~ f cos(v^v':^)cos(V2v3^)sin(26lp+); (56) 



efeFl^+fcF2tXfcF14, + fcF24, ~^ H e. 



(57) 



- fsin(%/2(p::^)cos(\/2(^^)sin(26'p+); (58) 
XfeFit+fcF2t<^fcFu+fcF24- + H-c. ^ (59) 
~ fcos(y2i^:^) sin(V2v3^) sin(26'p+) . (60) 

Here F = ?yi-t-?7i^772t'?2^- Note that we have only listed ob- 
servables containing siii(20p+). The other independent spin-0 

objects e^r and jri^' contain cos(26'p_|-) and vanish because of 
the pinning condition Eq. ( |33] |. 

Having discussed S*^ -conserving observables, we can sim- 
ilarly consider S'^-raising observables. We will call objects 
corresponding to SS'' — 1 or 2 as "spin-1" or "spin-2" respec- 
tively. We have spin-1 bilinears. 



t+P'kFb 



(61) 



Generically, these all carry different momenta. We can readily 
write bosonized expressions. For reference, we give the main 
ones that contain oppositely moving fields: 



(62) 



where we used convention a,b = +/— for band 1 or 2. We 
can generally argue that at Q ^ 0, tt, objects Sq that trans- 
form like Fourier modes of the 5'+(.t) operator cover, up to 
complex phases, all distinct spin-1 observables. In the present 
SBM system, we do not find any interesting spin-1 observ- 
ables at Q = 0, TT. 

Since we will encounter phases where 5"+ is gapped, we 
also need to consider SS" = 2 observables, i.e., some kind of 
"magnon pair" creation operators. Because of the hard spin 
condition, we define them on bonds [x, x + n], 

P+'(")(a;) EE S+{x)S+ix + n) (63) 
~ fl{x)fl{x + n)f^{x + n)f^ix) . (64) 

The last line can be expanded in terms of the continuum fields 
and organized as follows. For a-species, a "pair" opera- 
tor ,fa{x),fa{x + n) contains zero momentum contributions 
/iaa/laa' « = 1 or 2; ±(fcFiQ + fcF2Q) momentum contri- 
butions fpiafp2a^ P = L/R; and ±(fcFia ~ fcF2a) Contribu- 
tions fpiaf-P2a- Multiplying the pair creation operator for 
t species and pair destruction operator for l species, we ob- 
tain contributions to V'^ carrying combinations of the above 
momenta. 

We can argue on general symmetry grounds that, up to com- 
plex phases, there is a single independent spin-2 object at 
Q 7^ 0, TT. On the other hand, at Q = tt there are two indepen- 
dent objects that transform differently under lattice inversion; 
they can be reahzed by •p+'("-°™") ^jjjj p+-(n-odd) j.ggpg^_ 



tively. At Q = 0, we consider only objects 'Pgl^o ' 
the same symmetry properties for any n. 



which have 
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In the present SBM problem, the main spin-2 observables 
occur precisely at Q = 0, tt, and we give bosonized expres- 
sions only for these. For Q = 0, there are four possible terms: 



g-i[2ip„+ + \/2(ai/3^ -bip^ '''(65) 



with independent a,b = +/— corresponding to bands 1 or 2. 
For Q = TT we find 



Q = TT 



sm 



(66) 



Because of the pinning condition on the 6'p-|_, only the 
^Q=^ °'^'^^ are non-zero, and we can use the nearest-neighbor 
magnon pair operator 7^+'(i) as the main representative. 



V. NEARBY PHASES OUT OF THE SBM IN THE FIELD 

We now consider what happens when either w^^ or ^^12 
from Eq. ( |25] l or both become relevant. 



A. Phases when is relevant 

Let us start with the case when the w^^ '^^'"m is relevant 
while w\r^ is irrelevant. The field (pZ^ is pinned, while fields 

(yS^ and remain gapless, so we have two gapless modes. 
There is no static order We summarize characteristic power 
law observables in Table U and discuss them in turn. 

First, all observables eg and XQ in Eqs. (I48ll-(l52]i con- 
structed out of the fields show power law. On the other 
hand, such observables constructed out of the fields that 
contain Ql^ become short-ranged once we pin the conjugate 

<^^; thus, only Q =- kp\^ + kp2^ can remain power law. There 

are two cases depending on the sign of wjj- 



^ ^ _ (2n-|-l)7r 

> L) : — — :z~f^ — ' n ^ L ^ 



2^2 



_ 2n7r 



(67) 



(69) 



Next, note that all spin-1 observables S'g become short- 
ranged since they all contain the wildly fluctuating field . 
Schematically, the individual become gapped because of 
their "pairing". On the other hand, spin-2 observables contain 
pairs of and can remain gapless. Explicitly, after pinning 
the Lp^ , we have for the dominant correlations at Q = and 

TT 



' Q=0 



(71) 
(72) 



Pinned ip^ : Common power-law order for either sign of 




^±(fcFl^+fcF24,)' 
X±(fcFli+*:F2i) 






v+ 


Distinct power law correlations 




w'l^ > : 




) 
) 





TABLE I: Summary of the main observables when lujj t^rm is rele- 
vant and pins Lp^ . Critical wavevectors Q for the magnon-pair cre- 
ation operator are obtained by combining any of g-f^ = {0, ±(fci?i-|- + 
fcF2t)} withanyofg; = {0, ± (fc^i; + fcF2;), ± (^fu - A;F2i)}, 
Q ~ + Ql', the most important ones are Q — and n. 



The gaplessness of the Lpcr+ is required since S'fgj is con- 
served and incommensurate with the lattice. We can map the 
spin system to hard-core bosons^^ and in the present case 
single boson excitations are gapped, while pair boson exci- 
tations are gapless and created by e^'^'^''+ .... In the "particle- 
hole" sector, we have strong "density" or "current" correla- 
tions, Eq. ( |68] | or (iTOl l. at wavelengths that can be related to 
typical separations between boson pairs, and such e*^''+ con- 
tribution is generally expected in a Luttinger liquid of pairs. 
Thus, the resulting state has spin-nematic power law corre- 
lations as well as density or current power law correlations. 
Which one is dominant depends on the scaling dimensions of 
gi2v5,T+ versus e'^''+ . The scaling dimensions would need to be 
calculated numerically since the <^cr+ and (p^ mix in general; 
we do not attempt such quantitative estimates here. 

Having discussed observables controlled by the gapless a+ 
part, let us finally mention that ySi^'' and Xtt directly detect the 
gapless Lp'^ field, cf. Eqs. (I53]l-(l60ll. In the phase discussed in 
this section they have the same power law decays. 

We have considered the case when only becomes rele- 
vant. The case when only ■W12 becomes relevant can be treated 
similarly by interchanging ^ and |. 



Let us now discuss the phases out of the SBM when both 
w{2 and 'W12 terms get relevant. Once the couplings flow to 
large values, both variables p^ and p^ will be pinned so as 
to minimize the energy. There are four possible situations de- 
pending on the signs of the wjj and ■ 

In all cases, we find that the translational symmetry is bro- 
ken by either a static order in Bt^ (corresponding to period-2 
valence bond solid) or x-k (corresponding to period-2 chiral- 
ity order). Coexisting with this, we have one gapless mode, 
namely the overall spin mode "(7+", which must remain gap- 
less as long as the magnetization density is incommensurate 
with the lattice. Similarly to the case with one relevant cou- 
pling, spin-1 observables are gapped. Spin-2 observables are 
gapless, with the dominant contributions 



I Q=0 ' Q= 



,-i2(^„ 



(73) 



(Note that the original wavevectors Q = and tt are not 



t 


i 
"^12 


Static Order 


Power-Law Correlations 


+ 


+ 






■P + 

{Q} 








X±(feFlc+fcF2„) 


' {Q> 


+ 






^±(fcFlt + '=-F2t)' 
X±(feFli+feF2i) 


' {Q} 




+ 




^±('!FU + '=F24.)' 
X±(fcFlt"^^F2t) 


V+ 
' {Q} 



TABLE IL Summary of the cases when both to^j ™d terms are 
relevant. For Wj2'"^i2 > we have period-2 VBS order, while for 
"'12^12 < we have period-2 chirality order. In all cases, coexisting 
with such static order, we have power law correlations in the spin-2 
(magnon pair) observable and in the specific e/x observables. 

distinguishable once we have the period-2 static orders.) To- 
gether with such spin-nematic observables, we also have spin- 
observables of the e- or x-type depending on the pinning 
details, with the wavevectors ±{kFia + kp2a) which satisfy 
kpit + kp-it = -{kpu + kF2i) - Ti". 

Below, we consider four different pinning situations in 
more details. The main features in each case are summarized 
in Table nil 

1. wl^ > 0, TO^a > 
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FIG. 5: Picture of the Valence Bond Solid order when gains 
an expectation value. Top: Id chain view. Bottom: the same in 
2-leg ladder view. Coexisting with the static order, we also have 
spin-nematic power law correlations and power law in either e or x 
channels (these properties are not depicted in any way). 

FIG. 6: Picture of the static period-2 order in spin chirality when Xtt 

(2) 

gains an expectation value. Since jTif ~ Xt, we have static stag- 
gered second-neighbor bond currents in the chain view (top figure). 
In the ladder view (bottom figure), we have oppositely oriented spin 
currents flowing on the two legs. Coexisting with the static order, 
we also have spin-nematic power law correlations and power laws in 
e/x channels (these properties are not depicted in any way). 



The pinning conditions for fields ip^ and ip^ are 

(2n + l)7r (2m+l)7r 



2^/2 

3(1) 



2V2 



In this case, obtains an expectation value while Xtt ~ 0. 
Thus we expect to see period-2 VBS order as illustrated in 
Fig.|5] We also have power law correlations in 



^fcFlt+^F2t ^ ^ — 'i^'F14,— ^F2 



(75) 



while XfcFic,+fcF2c, = 0- Note that because of the relation 
Eq. {5^ [in the sense that i{ekj,^^+ki,2r<^kpn+kp2i~^-(^-) has 
the same symmetry properties as yBl^"*], once the system devel- 
ops static order in Bi^\ the ekpi-f+kp2r e-fcFu-feF2i = 
^kFij^+kF2i ^'^ longer independent. Appropriately, the 
wavevectors kpi^ + kp2^ and —kpii — kp2i differ by tt and 
also become connected. 



2. < 0, < 



Here, the pinning conditions are 

_ 2mT _ 2rmr 
(/^^ = fl=^r^i n,meZ. (76) 



2V2' 



Again, sl^^ obtains an expectation value while Xtt = 0. How- 
ever, here we have power law correlations in 

(77) 



while £kFia+kF2a = 0' Similar to the discussion in the 
preceding case and using relation Eq. dSSl l. XkFi-t+kF2t ^"'^ 
X-A;Fu-fcF24 are not independent observables in the presence 
of the static order in Bn 



(1) 



3. wl^ > 0, wfj < 

In this case, the pinning conditions are 

(2?! + l)7r _ 2m7r 
= — — ' = n,meZ. (78) 



2V2' 



In this phase, Xtt obtains an expectation value while Bi^^ = 0. 
Thus we expect to see period-2 chirality order as illustrated in 
Fig.|6l We also have power law correlations in 

£kpi^ + kp2-t ^ X-kpi^-kp2i ~ ^ ' (79) 

while XfcFit+feF2T = e-fcFu-fcF24, = 0. By using Eq. (|57]), 
we can understand the equivalence of the two observables 
ekpi-f+kp2r and x-kpn-kp2i once there is the static order in 

Xtt- 



4. < 0, luta > 

In this case, the pinning conditions are 

2?i7r _ (2m+l)7r 



2V2' 



2^/2 



(80) 
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Xtt obtains an expectation value while B)t 
power law correlations in 



(1) 



0. We also have 



(81) 



while ekp„+kf 



X~kFi^ — kF2i 



0. The two observ- 



ables XfeFiT+fcf2T and e_fc^ij^_fc^,^ 



become related because 
of Eq. (|59] | and the static order in Xtt- 

This completes our discussion of the phases out of the 
SBM. We cannot tell which of the different cases are more 
likely in particular microscopic models. Also, the power law 
correlation exponents depend on the unknown Luttinger pa- 
rameter of the "(T+" field, and we cannot tell whether 
spin-2 or spin-0 observables dominate (their scaling dimen- 
sions are l/ga+ and ga+f'i respectively). However, we have 
developed a qualitative understanding of the phases and ob- 
servables needed to identify them, which we hope will be use- 
ful in numerical studies of models reahzing the SBM phase. 



VI. DISCUSSION 

In this paper, we studied instabilities of the 2-leg SBM un- 
der the Zeeman magnetic field. The instabilities are driven by 
the interactions, Eq. dH), and we analyzed possible out- 
comes using Bosonization. In all cases, we found a gap to 
spin-1 excitations, while spin-nematic (two-magnon) correla- 
tions are power law. Loosely speaking, this appears because 
of some pairing of spinons, while the precise characterization 
is obtained by analyzing all physical observables. 

Here we want to discuss consequences if such spinon pair- 
ing were to occur in a 2D spin liquid under the Zeeman field. 
At present, we do not have any energetics justification under 
which circumstances this may happen and whether this ap- 
plies to the candidate spin liquid materials. However, the re- 
sulting states are quite interesting on their own and perhaps 
such phases may occur in some other 2D systems (several 
papers^^i^ considered mechanisms for spinon pairing in zero 
field). 

First of all, the analog of the stable SBM phase in Sec. |III] 
has gapless Fermi surfaces for both f and J, spinon species, 
with somewhat different kp'f and kp^. In the organic k- 
(ET)2Cu2(CN)3 and EtMe3Sb[Pd(dmit)2]2 materials, we es- 
timate (7i-|- — ni) I (n^ + n^) < 0.02 under laboratory fields, 
so the difference between the two Fermi surfaces is small. In 
mean field, the spin correlations are 



(5+(r)5-(0))^ 



cos[(kFt + k^i) • r + f ] 



cos[(kj^^ - k;^^) • r] 



(82) 



(83) 



l + cos[2k. -r+f] ^^^^^ 
z — ^ W\<* 



while gauge fluctuations are expected to enhance the k 
ki?^ and 2kj?ct partSf22, similarly to the ladder case^i^ 

Next, we want to discuss the analog of the situation in 
Sec. IV Al where there is pairing in one spinon species (say. 



/^) and no pairing in the other species. Note that the pairing 
must be odd- wave since it is within one fermion type. We will 
not consider any energetics selection of the pairing and just 
mention possibilities like p-wave {p^^ + i-py) or /-wave that 
can be nicely placed on the triangular lattice. 

The properties of the resulting phase are as follows. The 
species are gapless with Fermi surface, so we expect metal- 
like specific heat C — 7T; note that this is the full result since 
the gauge field is Higgsed out by the pairing. We also 
expect constant spin susceptibility at T since both 
and systems are compressible, the former due to the pair- 
condensate and the latter by virtue of finite density of states at 
the Fermi level. Because of the Fermi surface, we expect 
(5^ (r) 5*^(0)) to show 2kp\^ oscillations with Xjr^ power law. 
On the other hand, (5+(r)5"(0)) will show either a full gap 
if the /|- pairing is fully gapped as in the case of j)^ + iVy pair- 
ing, or a pseudogap if the pairing has gapless parts as in 
the case of /-wave pairing. Note that this does not contradict 
the finite susceptibility since the /^-pair condensate can read- 
ily accommodate AiV^ = ±2 changes. Related to this, spin- 
nematic correlations are gapless and show 1 /r^ power law at 
zero wavevector (in the mean field calculation). Interestingly, 
the gap or pseudo-gap to spin-1 operators would have conse- 
quences for NMR experiments done with ^H or ^■^C that are 
both spin-i nuclei and relax only by spin-1 excitations. From 
such measurements, this phase might appear gapped, but it 
actually has a gapless Fermi surface of one species. (In the 
context of ID models exhibiting spin-nematic phases, conse- 
quences for the NMR relaxation rate were discussed in detail 
e.g. in Ref.lM) 

Finally, let us consider the analog of the situation in 
Sec. IV B| where both /^ and /^ become paired, with possibly 
different pairing A]^^.,, A^^,. In this case, 5*^ and 5*+ correla- 
tions are both gapped (or pseudo-gapped), while spin-nematic 
correlation shows long-range order. Specifically, in the mean 
field. 



(5+(r)5+(r'))„,f = A:;,A^ 



(85) 



Note that this nematic order resides on the bonds of the lattice 
and details depend on the A^ and A^. For example, if we 
take A^ and A^ to have the same pattern, this will give ferro- 
nematic state. Curiously, if we take A^ ^ Vx^ ipy and A"^ ^ 
Px — iPy, we get q = antiferromagnetic nematic order on the 
Kagome lattice formed by the bonds of the triangular lattice. 
We emphasize that we have not discussed any energetics that 
may be selecting among such states. Whether something like 
this can appear in realistic models on the triangular lattice is 
an interesting open question. 
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